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Graphene research is currently one of the largest fields in condensed matter. Due to its unusual 
electronic spectrum with Dirac-like quasiparticles, and the fact that it is a unique example of a 
metallic membrane, graphene has properties that have no match in standard solid state textbooks. 
In these lecture notes, I discuss some of these properties that are not covered in detail in recent 
reviewfl^. We study the particular aspects of the physics/chemistry of carbon that influence the 
properties of graphene; the basic features of graphene's band structure including the tt and cr bands; 
the phonon spectra in free floating graphene; the effects of a substrate on the structural properties 
of graphene; and the effect of deformations in the propagation of electrons. The objective of these 
notes is not to provide an unabridged theoretical description of graphene but to point out some of 
the peculiar aspects of this material. 

PACS numbers: 81.05.ue, 63.22.Rc, 61.48.Gh, 73.22. Pr 



I. INTRODUCTION 

Graphene, a two-dimensional allotrope of carbon, has created a immense interest in the condensed matter commu- 
nity and in the media since its isolation in 2Q0^. On the one hand, graphene has unique properties that derive from 
its honeycomb-like lattice structure such as the Dirac-like spectrum (that mimics effects of matter under extreme 
conditions), its low dimensionality (that leads to enhanced quantum and thermal fluctuations), and its membrane- 
like nature (that mixes aspects of soft and hard condensed matter). On the other hand, because of the strength 
and specificity of its covalent bonds, graphene is one of the strongest materials in nature (albeit one of the softest), 
with literally none extrinsic substitutional impurities, leading to the highest electronic mobilities among metals and 
semiconductors^. Therefore, graphene is being considered for a plethora of applications that range from conducting 
paints, and flexible displays, to high speed electronics. In fact, it can be said that perhaps, not since the invention 
of the transistor out of germanium in the 1950's, a material has had this kind of impact in the solid state litera- 
ture. However, unlike ordinary semiconductors such as germanium, gallium- arsenide, and silicon, graphene's unusual 
properties have to be understood before it can really have an impact in technological applications. 

Any material has a hierarchy of energy scales that range from the atomic physics (~ 10 eV), to many-body effects 
IQ^^ eV). To understand the behavior of a material one needs to understand how these different energy scales 
affect its macroscopic properties. While structural properties such as strength against strain, shear and bending may 
depend on the covalent bonds formed by the atoms, magnetism and superconductivity are governed by the particular 
way electrons interact with each other through Coulomb forces. Furthermore, while the properties of metals and 
semiconductors depend on the physics close to the Fermi energy (a direct consequence of Pauli's exclusion principle), 
the nature of the vibrational spectrum depends on the particular way ions interact among themselves and how the 
electrons screen these interactions. 

One of the great accomplishments of the application of quantum mechanics to the theory of metals is the under- 
standing that while different materials can be structurally very different from each other, its long wavelength and low 
energy physics is essentially identical and depend on very few parameters. This so-called renormalization towards 
the Fermi energjE is one of the greatest theoretical accomplishments of the twentieth century and is the basis of 
Landau's theory of the Fermi liquicPl. In systems where the low energy physics is described by a transformation of 
coordinates with Galilean invariance the most significant parameter is the "effective" mass of the carriers that acts 
as to generate a scale from which is possible to compute most of the important physical quantities such as specific 
heat, magnetic susceptibility, electronic compressibility, and so on. The most basic difference between graphene and 
other materials is that its low energy physics is not Galilean invariant, but instead Lorentz invariant, just like systems 
in particle and high-energy physics with Dirac particles as elementary excitations®. In this case, the renormalization 
towards the Fermi energy is different from other material^ because, in the absence of a "mass" (which in graphene 
means the absence of a gap in the electronic spectrum), all physical quantities depend on a characteristic "effective" 
velocity that plays an analogous role as the speed of light plays in relativistic quantum mechanic^. However, unlike 
true relativistic fermionic system^, the Dirac quasiparticles in graphene still propagate with a velocity that is much 
smaller than the speed of light, the speed that Coulomb interactions propagate. Therefore, the Coulomb field can be 
considered instantaneous in first approximation, making of the electrodynamics of graphene electrons a mix between 
a relativistic and a non-relativist problem. Clearly this unusual situation requires a re-evaluation of the Fermi liquid 
theory for this material. 
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In these notes, however, we are going to focus on the more basic aspects of graphene's properties. In Section [TT] 
we are going to discuss the s — p hybridization theory and how it leads to the basic energy scales of the graphene 
problem. We are going to show that the problem of hybridization is controlled by the angle between the s — p 
hybridized orbitals. Using the Slater-Koster theory, we compute the main matrix elements of the problems in terms 



of this angle. In Section III we move from molecular orbitals to the crystal and discuss the simplest tight binding 
Hamiltonian that describes the full band structure, that is, that includes both tt and a bands. We show that even 
this simple band-structure reproduces quite well the results of more sophisticated ab initio methods. This particular 
description becomes particularly good close to the Fermi energy where the Dirac particles emerge naturally. Phonons 



in free floating graphene are discussed in Section IV We show how the flexural modes result from the bending energy 
of a soft membrane and how those modes can be quantized. We show that as a result of the presence of flexural modes 
the linear phonon theory predicts an instability of the graphene sheet towards crumpling. The effect of a substrate 
on the flexural mode spectrum is discussed in Section |V] We show that the presence of a substrate, that breaks 
rotational and translational symmetry, allows for new terms in the phonon Hamiltonian that change considerably the 
energy dispersion of the flexural modes. We also show that within the linear theory, graphene follows the substrate 
in a smooth way with the characteristic length scales that are dependent on the details of the interaction with the 
substrate. When graphene is deformed in some way, either by bending or strain, the electronic motion is affected 



directly. In Section VI we show that at long wavelengths deformations lead to new terms in the Dirac equation that 
couple to the electrons as vector and scalar potentials. There are cases, however, where graphene is not slightly 
deformed but strongly deformed in which case the Dirac theory has to be completely reconsidered. In these notes 
we discuss briefly the case of a graphene scroll which results from the competition between the bending energy (that 
favors flat graphene) and the van der Waals interaction of graphene with itself (that wants it to have maximum area 
overlap) . This is only one example where structural deformations can have a strong effect on many of the electronic 
properties of this material. Our conclusions are given in Section [VII[ 

II. THE CHEMISTRY 

The electronic configuration of atomic carbon is Is^ 2s^ 2p^. In a solid, however, carbon forms s — p hybridized 
orbitals. The Is electrons form a deep valence band and essentially all properties of carbon-based materials can be 
described in terms of the 2s and 2px, 2py and 2pz orbitals that can be written a,ff^. 

{r\s) = i?,(r) X 1, 

(r|p,) = Rp{r) X VS sin 6 cos 6 , 

(r|p.y) ~ i?p(r) X \/3 sin 6* sin ^ , 

{r\p,) = Rp{r) X 73 cos 0, (1) 

where Rs{r) = (2 — r/ao)e~^^^'^°'°\ and Rp{r) = (r/ao)e~''/^^°''^ are the radial wave-functions. One particular way to 
parametrize a hybridized s — p state is given by Pauling^-^: 

|0) - A\s) + ^l-A^\pz) , (2) 

where A is a parameter that describes the degree of hybridization between s and p states. This basic orbital is shown 
in Figjl] The energy associated with this orbital can be obtained from the hydrogen atom spectrum: 

Eo^e^^{0\H^\0)^A^Es + {l-A^)Ep (3) 

where iJo is the hydrogen atom Hamiltonian where Es ~ —19.38 eV is the energy of the 2s-state and Ep ss —11.07 
eV is the energy of the 2p-state. 

All the other 3 orthogonal orbitals can be constructed starting from ([2|. A particularly simple parametrization is 
the followin^i2l 

|1) = - ^')/3k) + v^b.) - {AlV^)\p.) , 
|2) = v/(i - A^)/2,\s) - yiTeb.) - VV^\Py) - {A/Vi)\p. 



|3) = V(l - A^)/5\s) - ^/^76b.) + Vij2\py) - {A/V3)\p.) . (4) 



Notice that A controls the angle between the z axis and these states. We can clearly see that the direction of largest 



FIG. 1: Absolute value of the wave-function for a hybridized s — p state, Eq.lpl), for A = 0.5. 



amplitude for one of these orbitals (say, (r|l)) is given by: 

^(r|l)(0 = O) = V2cose,^ + Asine,^^Q, 

= -arctan(%/2/A) . (5) 

Hence, for A = the hybridized state |1) is perpendicular to the other orbitals that remain in the x — y plane. This 
is the so-called sp^ hybridization (see Figj2|. For A = 1/2 the orbitals have tetragonal structure making an angle of 
109.47 degrees with the z axis. This is the sp^ hybridization (see FigS. 




FIG. 2: sp^ (left) and sp^ (right) hybridized orbitals. 



In free space, the states |1), |2) and |3) are clearly degenerate and their energy is given by: 

E^^E2 = E:i = (l|i/o|l> - - + ^"IT"^^' 

The energy of these states is shown in Fig. [s] Notice that in the sp^ case {A = 1/2) all orbitals are degenerated while 
in the sp^ case the orbitals are separated by an energy of approximately 2.77 eV. 

The presence of another carbon atom induces a hybridization between the different orbitals. This hybridization 
depends on the distance ^ between the atoms and also on the orientation of the orbitals relative to each other. The 
distance dependence is usually well described by an exponential behavior: 

K(^) « (7) 

where — d\Ti{Va) / dt where a labels the different orientations of the orbitals. In terms of orientation, there are four 
different types of elementary hybridization between different orbitals (shown in Fig|4|: Vssa (~ —5 eV for I = 1.42 
A); T4p<T (« +5.4 eV for ^ = 1.42 A); (w -h8.4 eV for € = 1.42 A); F^p^ (w -2.4 eV for £ = 1.42 A^. 
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FIG. 3: Energy of the hybridize s — p states, given by S and ([6| as a function of the hybridization parameter 
corresponds to the sp^ and A = 1/2 corresponds to the sp^^ configuration. 
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FIG. 4: Basic hybridization energies for s — p bonds: (a) Vsso-; (b) Vapa; (c) Vppa; (d) Vp, 



Any hybridization energy can be obtained from those basic hybridizations shown in Fig. [4] as linear combinations. 
For instance, let us consider the intra-atomic hybridization between between |2) and |3). This is given by the matrix 
element: 

l^„tra - {2\Ho\3) = (2| (J^^E^l-^) - ^\P.) + ^\Py) 



(Es - Ep) 



(8) 



We can also compute inter-atomic hybridization energies such as the hybridization between two |2) states oriented as 
in Fig. [5] 



spa ' 



(9) 



III. THE CRYSTAL AND BAND STRUCTURE 



In this section we focus on the flat graphene {sp^ only) lattice. The geometry of the lattice is shown in Figj?] (a) 
and it can be easily seen that the honeycomb lattice can be considered a crystal lattice with two atoms per unit cell. 
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FIG. 5: Calculation of the hybridization of two s — p orbitals in terms of the basic hybridization energies shown in Fig|4] 




FIG. 6: Honeycomb lattice and electronic orbitals 



These two sublattices, A and B, form two interpenetrating triangular lattices. The lattice vectors are: 



ai 
a2 



(3, VS) , 
(3,-^), 



(10) 



and the nearest neigh boor vectors are: 



Si = a(l,0) 
a 



S2 



(11) 



where a = 1.42 A is the carbon-carbon distance. The first Brillouin zone is shown in Fig. [rjjb) and the reciprocal 
lattice is spanned by the vectors: 



da 



27r 



(1,-V3). 



(12) 
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FIG. 7: (a) Lattice structure; (b) First Brillouin zone 



The K point, located at 



3a V V3 



(13) 



is of particular importance here. 

Let us label the state at each carbon atom as: \n, a, i) where n = 1, N labels the unit cell, a = A,B labels the 
sub-lattice (see Fig[7][a)), and i = 0, 1, 2, 3 labels the states defined in ([2| and Q. In the simplest approximation we 
can consider a tight binding description where the electrons hop between different orbitals (intraband) in the same 
atom with characteristic energy Vlntra, between the two different sublattices (interband) along the planar orbitals with 
energy V^, and between two different sublattices with orbitals perpendicular to the plane with energy V^p^. In this 
case the Hamiltonian can be written as: 



intra i'^nM,i'^"-,"-d ~^ h.C.) 



^n.A.i^'m.B,! + h.C.) , 



(14) 



(n,m) ,i^0 



The first step for the solution of this problem is to use the fact that the system is invariant under discrete translations 
and to Fourier transform the operators: 



k 



(15) 



which leads to: 



H = Yy^Yl 4,a,0Ck.a,0 + e<T XI 4,a,iCk,a.» + V^intra Yl (4,a,»Ck,aj + h.C.) 

+ l^ppTT 7k(cl,_^_oCk,B,o + h.C.) + K \e'^'^'ci,A,iCKB,i + e"^'^^ ,^ .^0,^3,2 + e''''*'c^ gCk.s^a + h.c. 



} (16) 



where 



7k = > e 

1=1,2,3 



E 



|7k| = V3 + 2cos(k-((5i-,52)) + 2cos(k-(5i-(53)) + 2cos(k-(^2-^3)). (17) 
In this case the Hamiltonian can be written as: 

^^ = E*k-[^]-*k, (18) 
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where Vfj^ — {c1^^a,0^ '^,B,Q^ ^L,A,1^ '^k,B,l^ ^k,A,2' '^,B,2' ^k,A,3^ ^^'^ 
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(19) 



It is immediately obvious that the Hamiltonian for the tt band completely decouples from the a band and can be 
treated separately. 



(a) 



Energy(eV) 





FIG. 8: Band structure of graphene (energy in ev). (a) 7r-bands; (b) cr-bands. 
Although (19) is a 8 X 8 matrix (a 2 x 2 block plus a 6 x 6 block), it can be diagonalized analytically to produce 8 



8 



energy bands that read: 



i?.,3,±(k) = + ^ - J (^)' + 1^,2 ± |yi,traT/.||7k| , (20) 



which represent two tt bands and six cr bands, respectively. The 7r-bands are shown in FigjSjJa) and the cr-bands are 
shown in Fig. [sj^b) . In Figjojwe compare the results of this simple tight-binding with more sophisticated calculationJSl, 
One can clearly see that the tight binding approach produces a fairly good description of the band structure, especially 
the TT bands, although more hopping parameters have to be introduced in order to reproduce the details of the a 
bands. The fact that this full band structure can be obtained analytically makes this particular parametrization of 
the bands rather attractive as a first step towards the electronic description of the graphene electrons. 




FIG. 9: Band structure of graphene (energy is given in Hartrees - 1 Hartree = 27.21 eV): left, result of the model discussed in 
the text; right, result of ref. [14| . 



As we said in Section |nj the carbon Is states are completely filled, leaving 4 electrons per carbon atom (8 per unit 
cell) to fill these 8 bands. Hence, the first 4 bands are fully occupied meaning that the Fermi energy crosses exactly 
at the "conical" points of the dispersion at the K and K' point of the Brillouin zone in Figj9j Notice that there are 6 
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of these points in the Brillouin zone with vectors given by: 



Qi = 4V(3%/3a)(0,l), 

Q2 = 4V(3\/3a)(y3/2,l/2), 

Q3 = 47r/(3\/3a)(\/3/2,-l/2) 



(21) 



and also at — Qi with i — 1,2,3. Close to these particular points, that is, close to the Fermi energy, we can find a 
simple expression for the electronic dispersion by expanding 7k as: 



7Qi+k ~ y(fe;r 



ik 



for k <^ Q. From (20) we see that the spectrum becomes: 



where 



E±{kx, ky) = ±VFk = ±VFJkl + fc2 , 



vp = ?,\VppTr\a/2, 



(22) 



(23) 



(24) 



is the Fermi velocity. Thus the electronic spectrum close to the Fermi energy has a conical form, mimicking the 
dispersion of a relativistic, massless, Dirac particle. Notice that in the first Brillouin zone each corner of the zone 
contains 1/3 of a Dirac cone but we do not need all the vectors in (21 1 to describe the problem. We can use reciprocal 
lattice vectors ( 12 ) in order to translate each piece of the cone to two corners located at iQi making 2 Dirac cones 
in the extended zone scheme. Hence, close to Qi we can rewrite the Hamiltonian as: 

kj^ 

k L 



y 







(25) 



where V'^.k = iQi+k and similarly for the B sublattice. By Fourier transforming the Hamiltonian back to real space 
we obtain the 2D Dirac Hamiltonian: 



Fn = y" rfV^'^(r)(wFCT- V-/i)*(r), 



(26) 



where n is the chemical potential measured away from the Dirac point, a = (ax, cry) are Pauli matrices, and 

*t(r) = (^^(r),^l,(r)). (27) 



IV. PHONONS IN FREE FLOATING GRAPHENE 



Since graphene is two-dimensional and has 2 atoms per unit cell there are 4 in-plane degrees of freedom that give rise 
to 2 acoustic, longitudinal (LA) and transverse (TA), and 2 optical, longitudinal (LO) and transverse (TO), phonon 
modes. However, graphene is embedded in a three-dimensional space and hence there are two extra degrees of freedom 
associated with out-of-plane motion: flexural acoustical (ZA) and optical (ZO) phonon modes. The displacements 
associated with these modes at long wavelengths is shown in FigjlOj In leading order in the displacements the in-plane 
and out-of-plane modes decouple and we can study them separately. Here we focus in the out-of-plane modes since 
those dominate the low energy physics of graphene. 

We use the Monge representation where the a point R in the graphene lattice is described by the vector R = (r, h{r)) 
where r — {x,y) is the 2D coordinate vector and h{r) is the height variable . The unit vector normal to the surface 
is given by: 

N = (-V/i + z)/v/l + (V/i)2, (28) 

where V — {dx,dy) is the 2D gradient operator, and z is the unit vector in the third direction. Distortions from the 
flat configuration cost energy because they rotate the sigma orbitals that want to be aligned in order to get maximum 
overlap of the wavefunctions that bind the atoms together. Let us consider the triangular lattice centered at sublattice 
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Transverse Optical (TO) 
*\ /* TOP VIEW 



Longitudinal Optical (LO) 




SIDE VIEW 

Flexural Acoustical (ZA) Flexural Optical (ZO) 

FIG. 10: Lattice displacements associated with some of the long-wavelength phonon modes in graphene. 



A, say at R^, that covers the entire graphene lattice. The normal at this point is N^. The energy cost for the bending 
of graphene can then be written as: 



(29) 



where is the lattice bending rigidity of graphene. Notice that (291 has the form of the energy of a classical 



ferromagnet where the normal vectors play the role of spins. Define Rj = + u where u is the next-to-nearest 
neighbor vector that connects two triangles and let us consider only length scales much larger than the lattice spacing, 
|u| w \/3a. In this, assuming that N is a smooth function over the lattice scale, we can write: 



and hence: 



N, «N, + (u. V)N, + -(u. V)2n, 



N, • N, « 1 + N, • [(u • V)N,] + -N, • [(u • V)2n,] 



where we used — 1. Again, assuming the distortions to be small and smooth, we can rewrite (28) as: 



N«z- V/i- -(V/i)2, 



(30) 



(31) 



(32) 



and hence. 



where a = x,y. Therefore, using (32 1 and (33) and the fact that z • V/i = 0, we find: 

N • [daN] « . 

Once again: 



and finally: 



N-[dl,N]^-{dbVh)-idaVh) 



(33) 



(34) 



(35) 



(36) 



Replacing (34) and (36) into (31) we find 



N, -Nj w 1 - - [(u- V)V/i] 



(37) 
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Replacing (371 in (29) we find: 



d^r[V^h{r)]^ , (38) 



where k (cx kl) is the bending rigidity. Notice that the energy (38 1 is invariant under translations along the z direction 
and rotations around any axis (the energy is invariant under h ^ h + Cia + €'2^ x 9 where Ci and C2 are c ont ains 
and 9 is the anti-clockwise oriented angle of rotation along a given axis). In momentum space we can rewrite (38 1 as: 

SUB^^j d\q%{ci)\\ (39) 

showing the elastic cost in energy for bending the graphene sheet behaves as when g — > 0. 

The quantum mechanics of bending can be easily obtained by quantization of the field h{r). Introducing a momen- 
tum operator P(q) that is canonically conjugated to /i(q): 

[Mq),P(q')]-*'5'(q-q'), (40) 

we can write the Hamiltonian for the bending modes as: 

^^|,.^|Phs^^^,(_q),(^)|^ (41) 

where a is graphene's 2D mass density. From the Heisenberg equations of motion for the operators it is trivial to find 
that /i(q) oscillates harmonically with a frequency given by: 

'^flcx(q)- (^)'^%^ (42) 

which is the dispersion of the acoustical fiexural mode. 

Notice that the existence of a mode dispersing like in two dimensions has some strong consequences. Consider, 
for instance, the mean square displacement of the height (w„ — 27rn//3 with /3 = \/T is the Matsubara frequency at 
temperature T): 



d^q 1_ 

(27r)2 (w„)2 + nq 



d^q 

(27r)2 2^q^ 

= / du ^ (43) 

where we had to introduce a ultraviolet cut-off (A ~ 1/a) and a infrared cut-off [l/L where L is the system size) 
because the integral is formally divergent in both limits. The ultraviolet divergence is not of concern since there is a 
physical cut-off which is the lattice spacing, the problem is the infrared divergence since it indicates an instability in 
the system. Let us consider two different limits. At low temperatures (T — >■ 0) we can approximate coth{u) w 1 and 
we get: 

(/.2(r^0))«-^ln(AL) (44) 

indicating that even at zero temperature the quantum fluctuations are logarithmically divergent. However, at high T 
we can approximate coth{u) sa l/it, the ultraviolet cut-off becomes irrelevant, and we find 

(/^'(T^^^))-^ (45) 
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and one finds a severe infrared divergence of the thermal fluctuations. Similar arguments can be made about the 
fluctuations of the normal N, namely, 

((5Nf)«((V/if)«^ln 

and hence at low temperatures we find ((^N)^(r 0)) « A^/ {Wtt^/k) and the problem is free of infrared divergences. 
At high temperatures we have (((5N)^(T oo)) « ln[AL]/(2/3K) and the fluctuations of the normal are logarithmically 
divergent. The interpretation of these results are straightforward: at zero temperature the graphene sheet should 
be rough but mostly flat, as one increases the temperature the thermal fluctuations become large and the graphene 
crumples (divergence of the normal indicates the formation of folds and creases). These results are modified by 
non-linear effects that we do not consider her^^^. More important than those, is the fact that in most graphene 
experiments, where the measurement of electric properties depends of physically constraining the samples, the basic 
symmetries of free floating graphene are explicitly broken. This is what we consider in the next Section. 



sinh(/3y;^A7(27r)) 
sinh(/3V^/(27ri2)) 



(46) 



V. CONSTRAINED GRAPHENE 



Since graphene is not floating in free space and since there can be impurities that hybridize with the carbon atoms 
changing the structural properties, it is important to consider symmetry breaking process that change the picture 
presented in the last Section. Consider, for instance, if graphene is under tension. In this case, the rotational symmetry 
along the x and y axis is broken and a new term is allowed in the energy, that reads: 

Ut = II d\ iVh{r)f , (47) 
where 7 plays the role of the surface tension. It is easy to see that the dispersion is modified to: 



io{q) = qJ-q^ + ^, (48) 
V cr <T 

indicating that the dispersion of the flexural modes becomes linear in g, as g — > 0, under tension. Notice that this is 
the case of a 3D solid, that is, the flexural mode becomes an ordinary acoustic phonon mode. This is what happens in 
graphite where the interaction between layers breaks the rotational symmetry of the graphene layers. It is easy to see 
that in this case the graphene fluctuations become suppressed: at low temperatures we find {h?{T — ^ 0)) is constant 
while (ft.^(T — ^ 00)) « rin(AL) is log divergent. Also, (((5N)2(T)) is well behaved at all temperatures. 

In the presence of a substrate, the translational symmetry along the z direction is broken and a new term is allowed: 

Us^lj d\{h{v) - s{v)f , (49) 

where v is the interaction strength with the substrate and s(r) is a reference height. Notice that in this case the 
flexural mode becomes gapped (by means of the h? term) and hence all fluctuations (quantum and thermal) are 
quenched. Therefore, the most generic energy for small and smooth height distortions of graphene is: 



^^2 



^ / d2r-|«[v2/^(r)]'+7[VMr)]'+t^(/^(r)-s(r))2} . (50) 



Let us consider the case of ( 50 1 in the presence of a substrate^fi. Minimization of the free energy with respect to h 
leads to the equation: 

KW^h-^iW'^h + vh^vs. (51) 
A particular solution of the non-homogeneous equation can be obtained by Fourier transform: 



where 



1/2 

1/4 
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are the length scales associated and curvature, respectively. Notice that (52 1 implies that the surface height more or 
less follows the substrate landscape, as expected. 

Let us consider the case of a random substrate where the probability of a substrate height between s and s + ds is 
given by: 



where A/" is a normalization factor: 



+CXD 



dsP(s) 



N 



(54) 



(55) 



where N is the number of surface cells, a is their lattice spacing, and sq is the average height variation. In this case, 
the height correlation function can be shown to be: 



s(r)s(r') = SQd{r — r') 



In momentum space (56 1 is written as: 



(56) 
(57) 



s(q)s(q')-(27r)2s4<5(q + q'). 
Using (57) we can immediately compute the height-height correlation function as a function of the in-plane distance: 



hir)h{0) 



d'q f d^q' 



2n 



dq 



where we have defined: 



(1 + £2q2 _^ ^4^4)2 2nif 



(27r)2 J (2^)2 (1 + ^2^2 + ^4g4)(l + £2(^,)2 + £4(^,)4) 

oo 



—^s{(\)s{q!) 



dq 



qJo{qr) 



(g2+p2)2(^2+p2)2 



(58) 



C = 2 



7 ' 
= 2 



(59) 



and the behavior of the integral depends on whether ^ is smaller or bigger than 1. For r S> A = l/Min(Re(P±)) one 
finds: 



(60) 



showing that the height fluctuations are short ranged and decay with a characteristic length scale given by A. The 
local height variations are given by: 



where for Q > 1 we have: 



and for C < 1- 



^-(C) 



(C^ - If'^ |-arctanf^^ 



(1 - C2)3/2 



1- v/W^ 



(61) 



(62) 



(63) 



the function w{C,) is shown in Fig. 11 Notice as C ^ 1, that is, when ^ ^t, the height variations are strongly 
suppressed and the system is essentially flat. This happens because in this limit the interaction with the substrate is 
strong, and the bending rigidity is large, compared with the tension. 
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FIG. 11: Amplitude of the local height variations in graphene supported in a random substrate, w{(). 



VI. DEFORMED GRAPHENE 



Graphene is a one atom thick and hence very soft. Just as any soft material, graphene can be easy deformed, 
especially in the out of the plane direction where the restitution forces vanish. These deformations couple directly to 
the electrons, making of graphene a unique example of a metallic membran^^JD, Lgt us consider what happens to the 
electrons when the graphene is deformed either by bending or strain. One effect is the change in the distance between 
the atoms, the other is the change in the overlap between the different orbitals. In both cases, the hopping energy 
between different carbon atoms is affected. Consider the case that of the nearest neighbor hopping that changes from 
Iq to tg + Sti{S) at site R,; in the direction of S. In this case we have to add a new term to Hamiltonian (16): 



5Ho = -;^45k'E'^^^w^'^'^"''^''-""'' 



where 



(64) 



(65) 



Notice that if we replace (64 1 into (26) the full Hamiltonian becomes: 



H = J (fr¥{r) a ■ [ivF^ + Aj - n ^'(r) , (66) 

where A — (Ax, Ay) and hence the changes in the nearest neighbor hopping couple like a gauge field. One should 
not worry about a possible broken time reversal symmetry here because there is another nonequivalent point in the 
Brillouin zone, at K', which is related to the K point by inversion symmetry. Hence, in the K' point the "magnetic 
field" is reversed and therefore there is no net time reversal symmetry breaking. 

Let us consider the problem of the hopping between next nearest neighbor sites. The Hamiltonian in this case is: 



Hi 



-tn 



{id) 



h.c. 



(67) 



which describes essentially hopping on a triangular lattice with lattice spacing \/3a with nearest neighbor vectors: 

vi = V3a(0,l), 



By Fourier transforming (67) we find 



V2 = V3a{V3/2,-l/2), 
vh = %/3a(-%/3/2,-l/2) 



(68) 



(69) 
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where 



6(k) = -t'oT.'^''"' 

V 

= 2<g |^cos(V3afcj,) + 4cos(3afc^/2) cos(\/3afcy/2) 



For k = Qi + q we find: 



e(Qi+q)«-3to 



Q^a^ 2 



(70) 



(71) 



for q <^ Qi. Notice that the first term ([71| leads to a shift in the chemical potential and the second term introduces 
a quadratic term in the dispersion. For q <C io/(3aio) '-^^^ disregard this term and consider only the chemical 
potential shift. Ab initio calculations estimate that tg O.lio providing a good range where this approximation is 
valitP. 

Consider again the deformations of the graphene surface. In this case in (71 1 has to be replaced by t' — + dt' 
and, in complete analogy with (64), we find: 



6Hi 



where 



(r) {ij\{r)ijA{r) + ^lir)ijB{i 



(72) 



(73) 



where we have used (111, (21), and ( |68[ ). In this case, changes in the next-to-nearest neighbor hopping couple to the 
Dirac fermions as a scalar potential. The final Dirac Hamiltonian has the form: 



H 



^(r) |(T • (^ivpy + Aj -1-1 + ^(r)! *(r 



(74) 



Given the local changes in the hopping energies the scalar and vector potentials can be readily computed through 
(65) and (73), respectively. 

This changes in hopping energies can be connected to changes in the structure of the lattice. Let us consider the 
case of in-plane distortions. In this case, the only change in the hopping energy is due to the change in the distance 
between the orbitals. Consider the case where the distance between sublattices changes by a distance S£ in the 
direction of S, in first order we have: 



dtK {dt/da)Si, 
(<5- V)u, 



where u{r) is the lattice displacement. Replacing the above expression in (65) we get: 

^^"^(r) « aa (u^^ - Uyy) , 

where a is a constant with dimensions of energy and we have used the standard definition of the strain tensor: 

1 / dui duj 



2 \dxj dxi 

An analogous calculation for the change in the next-to-nearest neighbor hopping energy leads to: 

4'("'(r) w g{u^x+uyy) . 



(75) 
(76) 



(77) 



(78) 



(79) 



Equations (77) and (79) relate the strain tensor to potentials that couple directly to the Dirac particles. 

Less trivial is the coupling to the out-of-plane modes since those involve rotations of the orbitals. In Fig{T2] we 
show a rotation of two orbitals by an angle 9. The rotation mixes tt and a states and for small angles the change in 
the hybridization energy is given by: 



t « to + ^ve"^ 



(80) 
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where 5V is the energy mixing between tt and a states. Notice that 9 = 211 /R where R is the curvature radius. In 
terms of the height variable we can write, in analogy with (371, that the change in the hopping amplitude due to 
bending in the direction u is given by: 



5t « 6V[{\i-\J)\Jhi 



(81) 



On the one hand, if u is a nearest neighbor vector 5 we get from (65 1 : 



{dlh + d^h) d^hdyh , 



(82) 



where the coupling constant Eat depends on microscopic details. On the other hand, if u is the next-to-nearest 
neighbor vector we find, in accordance with (73): 



(83) 



where Eaa is an energy scale associated with the mixing between orbitals. The main conclusion is therefore that 
for smooth distortions of graphene due to strain or bending the Dirac particles are subject to scalar and vector 
potentials leading to an "electrodynamics" that is purely geometrical (there is no electric charge associated with the 
"electric" and "magnetic" fields created by structural deformations). This structural "electrodynamics" has strong 
consequences for the electronic motion in graphene leading to many unusual effects that cannot be found in ordinary 
materials. In particular, one can manipulate the electrons by "constructing" appropriate deformations of the lattice 
that mimic electri c and magnetic fields. This is the so-called strain engineering and is a field of research that is still 
in its infancjU^I^. 




cos(e/2) + CiX3 sin(e/2) 



cos(e/2) + ( 



: siri(e/2) 



= Xp„coste/2)-Vpp„siri(e/2) 
^ V -(V +V )^ 



FIG. 12: Rotation of pz orbitals. 



A. A non-trivial example: The ScrolPSl 

So far we have discussed small distortions of the graphene sheet, but as a soft material graphene can be bent by 
large angles making completely new structures. The scroll is an example of the softness of graphene that shows up as 
an interplay between the van der Waals energy, Eyj^i^, that makes graphene stick to itself and the bending energy, Eg, 
the energy cost to introduce curvature in the graphene sheet. Let us show first that the scroll is a stable configuration 
of a free standing graphene sheet. 

We assume here that when graphene forms a scroll the equilibrium distance between folded regions is the same 
that one would find between the planes in graphite, that is, d = 0.34 nm. In this case, the scroll is described by the 
so-called Archemedian spiral: 



r{e) = ad, 



(84) 



where a = d/{2Tr). We assume that the scroll has length L along the scroll axis and consider the two main energy 
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FIG. 13: Graphene scroll in the form of an Archemedian spiral. 



scales in this problem: 



dW 



L 

Eb 

17 



Si 

-1 I ds{e) 

/eo+27r 



ds{e) 



1 



(85) 
(86) 



where 7 « 2.5 ev/nm is the van der Waals coupUn; 
the Archemedian spiral is given by: 



R{e) 



I K. 1 eV is the bending energji^. The radius of curvature of 
"VT , (87) 



and 



is the infinitesimal arc-length of the spiral. The evaluation of the integrals is straightforward: 



vdW 



7a 



In 



l(9 + 86|2) (6)0 + 27r)(9 + 8(6*0 + 271)2) 



{e„ + 2n) + ^l + {2n + eor {I + 91)^^ (1 + ((?o + 27r)2)3/2 



Eb 
L 



KG 
Y 



In 



'9i + ^i + ei 



(89) 
(90) 



9o + yrT 

The angle 9q is determined by the condition that the force on the graphene sheet vanishes at the edge of the scroll: 



1 dE 



(91) 



and using the values of the parameters we find 6q ~ 5.85. The total energy as a function of 9i is given by (89) and 



(90) and shown in Fig 14 Notice that the scroll becomes stable for angles of rotation bigger than sa 17.23 which is 
equivalent to approximately 3 turns of the scroll. This result is in agreement with more involved calculation^^. 

Consider what happens when a uniform magnetic field is applied perpendicular to the graphene plane. Notice that 
in the region of the scroll the component of the field perpendicular to the graphene surface is changing periodically. 
In fact, the normal component to graphene sheet at the scroll can be written as: 



N(0) 



1 



.(- sin(6l) - 6'cos(6l), cos(0) - 6'sin(6')) , 



and hence the component of the field perpendicular to the graphene is: 

B 



B^{9) = B • N(e') 



Vi + 02 



(cos(6') -6'sin(6i)), 



(92) 



(93) 



18 



6 
4 
2 

E(eV) 


-2 
-4 



12 14 16 18 20 

FIG. 14: Energy of a scroll as a function of angle of rotation. 

which is shown in Fig. |15| Notice that from the point of view of the electrons, which are constrained to Hve in a 
two dimensional universe, a uniform magnetic field in three dimensions, becomes an oscillating magnetic field in the 
presence of the scroll. An oscillating magnetic field has unusual effects on the electronic motion because the Lorentz 
force changes sign in the region where the field changes sign producing a "snake-like" motion around the reg ions of 
zero field. These snake states are chiral one-dimensional correlated systems with rather unusual properties^Sj. While 
it is very hard to produce a magnetic field that oscillates in a short length scale, it is very hard to roll graphene on a 
short length scale. Therefore, snake states should be easy to observe in rolled graphene sheets. 

1.0 
0.5 

B^fB 0.0 
-0.5 
-1.0 

10 20 30 40 50 60 

e 

FIG. 15: Magnetic field perpendicular to the graphene close to the scroll. 



VII. CONCLUSIONS 

Graphene is a unusual system that shares properties of soft and hard matter and mimics problems in particle 
physics. These properties have their origin on the nature of its chemical bonds, the s — p hybridization, and the low 
dimensionality. All these ingredients create a new framework for theoretical exploration which is still very much in 
its beginning. The possibility of studying electronics of deformed surfaces is quite intriguing because of its exotic 
"electrodynamics". The field of electronic membranes was born out of the discovery of graphene in 2004 and much 
has still to be understood. 

In these notes 1 covered some very basic aspects of graphene chemistry and physics. We have studied the basic 
Hamiltonian that describes the overall band structure of graphene. We have seen that flexural modes, nonexistent 
in three dimensional solids, are fundamental for the understanding of the structural stability of the material. We 
have seen that substrates can change considerably the flexural modes and hence control the height fluctuations in 
supported samples. Finally, we have studied the various ways that graphene deformations, either by strain or bending, 
can modify the electron propagation in this material. In particular, new structures that can be created out of flat 
graphene, such as the graphene scroll, can have exotic properties, such as snake states, in the presence of applied 
uniform magnetic fields. In some way, the present notes complete some of the material already published in ref. [T] 
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Nevertheless, the whole subject of graphene electronics and structure is much bigger than that and there is so much 
more to be understood. 

It is quite obvious that a material that is structurally robust, still flexible, and extremely clean, has enormous 
potential for technological applications. By understanding these unusual properties one can harvest new functionalities 
that did not exist before. Progress in material science unavoidably leads to progress. Graphene is considered today 
one of the most promising candidates for a new era of a carbon-based technology, possibly supplanting the current 
silicon-based one. However, only progress in the understanding the basic properties of these systems can actually 
drive the technological progress. 
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